Axisymmetric and triaxial MOND density-potential pairs

We present a simple method, based on the deformation of spherically symmetric potentials, to construct explicit axisymmetric and triaxial MOND densitypotential pairs. General guidelines to the choice of suitable deformations, so that the resulting density distribution is nowhere negative, are presented. This flexible method offers for the first time the possibility to study the MOND gravitational field for sufficiently general and realistic density distributions without resorting to sophisticated numerical codes. The technique is illustrated by constructing the MOND density-potential pair for a triaxial galaxy model that, in the absence of deformation, reduces to the Hernquist sphere. Such analytical solutions are also relevant to test and validate numerical codes. Here we present a new numerical potential solver designed to solve the MOND field equation for arbitrary density distributions: the code is tested with excellent results against the analytic MOND triaxial Hernquist model and the MOND razor-thin Kuzmin disk, and a simple application is finally presented. Subject headings: gravitation — stellar dynamics — galaxies: structure — methods: analytical — methods: numerical